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ABSTRACT 


This thesis investigated the current patterns around San Clemente Island, 
California. The data were from a mooring located at Eel Point, during the period July 10, 
2006 to July 23, 2007. The currents were dominated by strong poleward flow along 
isobaths. Semidiurnal and diurnal tides dominated the kinetic energy and rotary spectra. 
Clockwise rotation dominated the rotary spectra with inertial peaks in both the near 
surface and near bottom depths. There were on average 2-3 current reversals every three 
months which appear to be unrelated to the wind stress curl. The stratification was 
described for the year, with one upwelling event occurring in late spring. Using wavelet 
analysis, it was shown that inertial/diurnal and semidiurnal energy were present 
intermittently throughout the water column. The strongest modes of energy were 
associated with the near bottom depths. Wind stress was poorly related to the currents. 

Oscillations with 2-3 day period were observed in the alongshore flow. These 
lasted for approximately 10 days in April 2007 and are referred to as the April event. 
This event was observed from mid-water to the bottom. The event had maximum 
velocities of 30 cm/s at mid-water and an upward phase speed of 0.12 cm/s. These 
oscillations were rectilinear at mid-water and became counterclockwise with depth. It 1s 
speculated that this event may have been a coastally trapped wave but, with these data it 
cannot be determined as such. 

SCI/SCORE is owned and operated by Navy and supports training for the largest 
concentration of naval forces inthe world. This thesis provides information on ocean 
currents and stratification on the continental shelf at SCI for SCORE and military 
participants to use with military exercises, regional climatologies, ocean engineering 


design studies, and marine mammal studies. 
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I. INTRODUCTION 


Current and future naval operations are focused on coastal/littoral regions around 
the world. An accurate description of current variability is important to the Navy for a 
number of reasons, such as: mine-hunting, acoustic propagation, search and rescue, 
dispersion, small boat, amphibious and special operations. The United States Navy has 
owned, operated and trained on San Clemente Island, California since 1934. The 
Southern California Offshore Range (SCORE) is a division of Naval Air Force Pacific 
and is currently assigned operation control for all operations in and around San Clemente 
Island. The SCORE provides tactical range training and testing services to U.S. Navy 
units of the Pacific Fleet. The purpose of this study was to analyze and describe the 
coastal current variability on the continental shelf near Eel Point, San Clemente Island 
(SCI). SCI located in the Southern California Bight (SCB) region, 75 miles northwest of 


San Diego, California. 


A year of data from an ongoing study was collected in order to observe the 
currents near Eel Point. A mooring (Fig. I-l1) was located in about 130 m on the 
continental shelf at 32° 55’ N, 118° 35’ W. The objective was to gain a better 
understanding of current variability and stratification along the shelf of SCI. Due to 
Navy training exercises there have been few studies conducted at this location. This 
study investigates the physical forcings which may have played a role in driving the 
currents. Little is known about the circulation patterns around SCI and what is driving 
them. There are several studies which attempt to correlate the wind stress and ocean 
currents within the SCB (Winant and Dorman, 1997; Hickey, 1998; Bray et al., 1999; 
Hickey et al., 2003). Winds are a major driving force of the currents. This study attempts 
to relate wind stress and wind stress curl to the currents. This study will first describe the 
energy associated with the raw currents by calculating the principal axis of rotation, 
spectral analysis (kinetic energy and rotary spectra), and tidal analysis. In order to 
investigate the seasonal variability, the data set was divided into seasonal 3 month 


segments to describe the vertically averaged currents, along and cross-shore components, 


| 


and stratification throughout the water column. A wavelet analysis was run using the 
baroclinic currents and wind stress (calculated from three NDBC buoy stations, discussed 
later). The curl of the wind stress was calculated and then compared to the alongshore 
currents in an attempt to explain the current reversals observed throughout the time- 
series. A simple slab model following Pollard and Millard (1970) was used to investigate 


the oscillations observed in the currents. 


A. LITERATURE REVIEW 


1. Large Scale Circulation 


The large-scale ocean circulation that occupies the Northern Pacific basin 1s 
known as the North Pacific Subtropical Gyre. It is has a dominant anticyclonic flow 
which is comprised of four currents systems: the North Pacific Current to the north, 
California Current to the east, North Equatorial Current to the south and the Kuroshio 
Current to the west. The dominant patterns of the atmospheric circulation that drive the 
California Current System (CCS) are the North Pacific High, the Aleutian Low and in 
summer the thermal low over the southwestern United States. SCORE 1s located with the 


California Current System. 
Zz California Current System 


The California Current System (CCS) is an eastern boundary current system that 
runs along the west coast of North America from the Strait of Juan de Fuca to the tip of 
Baja California. The CCS is comprised of the equatorward-flowing California Current 
(CC); the poleward-flowing Davidson Current (DC), the poleward-flowing California 
Undercurrent (CU), and the Southern California Eddy (SCE). Seasonal variations in 
alongshore wind stress, coastline irregularities, bottom topography, and stratification, 
have been shown to be the mechanisms responsible for the observed large-scale structure 


within the CCS (Hickey, 1998). 


The CC is a slow (~2-10 cm/s), broad (~1000 km), equatorward year-round 
surface current. It carries colder, fresher subarctic water equatorward along the west 


2 


coast (Hickey, 1998). The current is strongest at the sea surface and generally extends 
over the upper 500m of the water column. The CC has a seasonal maximum in summer 
to early fall (Hickey, 1998) and 1s linked with seasonal jets and eddies. South of Point 
Conception, Lynn and Simpson (1987) use dynamic height fields to show the CC turns 
eastward, flowing into the Southern California Bight at 32° N. The eastward flow creates 
a cyclonic gyre (Southern California Eddy) within the bight that is strongest in late 


summer and fall. 


The CU 1s a narrow (~10-40 km) and relatively weak (~2-10 cm/s) current that 
flows poleward over the continental slope from Baja California to at least Vancouver 
Island (Hickey, 1998). The current is strongest at depths of ~100-300 m from the 
surface. It transports warmer and saltier subtropical waters poleward. Hickey (1979) 
concluded that the location, strength and core depth have shown considerable seasonal 
variability which can be related to the seasonal variability in wind stress and wind stress 
curl. It has seasonal maxima in late spring and summer with peak speeds of ~30-50 cm/s 


(Hickey, 1998). 


The DC flows poleward 1n the fall and winter from (~35° N) Point Conception to 
Vancouver Island (~50° N) (Hickey, 1998). It is a relatively broad (~200 km) and, weak 
(~5 cm/s) surface current. The DC seasonal maximum is in the winter. In Huyer et al. 
(1989), the reversal of winds from northwesterly in the summer to southeasterly in winter 
caused downwelling at the coast. This may be one of the forcing mechanism of the 


poleward flowing surface current known as the Davidson Current (Huyer et al., 1989). 
3. Circulation in the Southern California Bight 


The Southern California Bight (SCB) has a concave curvature between the 
headland at Point Conception and the U.S.-Mexican border (Lynn and Simpson, 1990). 
Due to the SCB’s complex bathymetry, it is considered to be a difficult region to 
determine seasonal circulation patterns. The bight is punctuated by a series of complex 
topographic features such as islands, basins and ridges. These exert strong control on 
circulation patterns at every depth and are governed by seasonal patterns in the CCS 
(Hickey, 1979, 1998, Lynn and Simpson, 1990; Bray et al., 1999, Hickey et al., 2003). 
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The CC flows equatorward offshore of the bight. A portion of the CC turns inshore and 
poleward, it then recirculates within the bight forming the Southern California 
Countercurrent (SCC). This circulation pattern has been referred to as the Southern 
California Eddy (SCE) (Figure I-2). The SCE is strongest in summer and early fall and 
weakest or even non-existent in spring (Hickey et al., 2003). The poleward flowing 
California Undercurrent (CU) lies beneath the sea surface (~100-200 m core depth) and is 


a dominant feature of circulation in the bight (Hickey et al., 2003). 


Spring and summer winds off Point Conception are strong due to the large 
pressure gradient between the North Pacific high and the thermal low over the southwest 
U.S. (Winant and Dorman, 1997). During the fall and winter, the atmosphere is 
dominated by traveling cyclones and anticyclones. Winds within the SCB (near the 
coast) are generally weak and directed onshore during the spring and summer. During 
the summer, the dominant wind direction 1s equatorward along the coast. This pattern 1s 
mostly influenced the by the North Pacific High. However in the winter it tends to 
reverse direction along northern California. Maximum equatorward wind stress occurs 
during the spring in southern California, and during late fall summer in northern 
California. Minimum equatorward stress occurs during late fall and winter (Hickey, 


1998). 


B. THESIS ORGANIZATION 


For this study, a one year period of ocean current observations from July 10, 2006 
to July 23, 2007 were used. The location of the of mooring site was near Eel Point which 
is located on the west coast of San Clemente Island (SCI), California. This thesis is the 
first to describe the current variability at this site. We will look first at the variability and 
stratification within the water column on the continental shelf (Eel Point mooring). Data 
collection will be discussed chapter IJ. Description of the currents and data analysis will 
be discussed chapter III. Chapter IV discusses the summary of results from the data 
analysis. Conclusions, future recommendations and Navy relevance are discussed in 


Chapter V. 
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Figure I-1.. | San Clemente Island, California and the position of Eel Point Mooring 
(red star) and NDBC stations 46025 (SMB-blue dot), 46047 (TB - green dot), 46086 
(SCB — orange dot). 
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Figure I-2._ | Major currents affecting the circulation of the SOCAL OPAREA. 
(Department of the Navy, 2005) 


I. DATA COLLECTION 


This chapter will describe data collection and processing for the moorings at Eel 
Point, San Clemente, California (Figure IIJ-1). The Eel Point mooring data used for this 
study consisted of two consecutive 6-month deployments SCORE 2 (SC2) and SCORE 3 
(SC3). This created a year long time-series. Data were collected by an upward-looking 
RD Instruments WorkHorse 300 kHz broadband Acoustic Doppler current profiler 
(ADCP) and three Sea-Bird Electronics SBE 37-IM MicroCATs. The Naval 
Postgraduate School (NPS) Oceanography Department and the Research Vessel Point 
Sur (R/V Robert Gordon Sproul) were responsible for deploying and (recovering) the 
moorings in July 2006 (January 2007). SC3 was recovered by R/V Point Sur on July 23, 
2007. 


A. EEL POINT MOORING (SC2 AND SC3) 


The Eel Point mooring was located on the west coast of SCI. Its location was 
chosen in consultation with NPS and SCORE so as to not interfere with naval training 
activities. Deployments of the mooring began in January 2006 with (SC1). For this 
study, we will use the data from deployments SC2 and SC3. SC2 was deployed on July 
10, 2006 at 32° 55.718’N, 118°36.503’W and recovered on January 18, 2007. The bottom 
depth at the mooring was 130m. A 300 kHz ADCP was deployed at a depth of 118 m on 
the mooring. Three MicroCATs were also deployed on the mooring at depths of 40 m, 
80 m, and 118 m. SC3 was deployed on January 18, 2007 at 32° 55.805’N and 118° 
35.676’W and recovered on July 23, 2007. The bottom depth at the mooring was 135 m. 
The ADCP was deployed at 120 m. The three MicroCATs on this deployment were at 45 
m, 85 m, and 123 m. Both deployments will be referred to as the Eel Point Mooring. 
Prior to the deployments, the ADCP and the MicroCATs were calibrated and tested by 
the NPS Oceanography Department. 


The instruments (ADCPs and MicroCATs) used to collect data at the Eel Point 
mooring were not set to sample at the same rate. The ADCP sampled every 15 minutes 
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and the MicroCATs sampled every 5 minutes. The RDI WorkHorse broadband ADCP 
uses an acoustic pulse and the Doppler shift principle to measure current velocities 
throughout the water column. The ADCP transmits an acoustic signal (a “ping’’) at a 
fixed frequency of 307.2 kHz. The transducers listen for the return echo from the 
backscatters in the water. The ADCP uses a four-beam transducer head with a 20 degree 
beam angle. In order to calculate the current velocities using a multiple beam solution, 
the ADCP assumes the currents are homogeneous throughout the horizontal water plane. 
Trigonometric and geometric relations between the beams are used to calculate the three- 
dimensional current velocity vectors. These are the u (north-south), v (east-west), and w 
(up-down) components. The fourth beam is a redundant beam used to calculate the error 
velocity and to the check the assumption of horizontal homogeneity. Error velocity is the 
difference between a velocity measured by one set of three beams versus a velocity 
measured by the other set of three beams during the same time frame (Simpson, 2001). 
Large differences could be associated with a non-homogenous environment or corrupt 


data caused by a large fish or obstruction blocking the transmission of the signal. 


The ADCP was setup with 4 m bin spacing and measures the average velocity 
over the depth range of each depth cell to produce a current velocity profile. Averaging 
over the range of each depth cell reduces the affects of spatial aliasing (Gordon, 1996). 
The profiles are generated by range-gating the echo signal (Gordon, 1996). Because 
shorter ranges return the echo signal faster than longer ranges, range gating was 
developed to break the received signal into successive segments for independent 
processing (Gordon, 1996). The ADCP was configured to burst sample. The burst 
sample consisted of 40 consecutive one-second pings during the initial 40 seconds of the 
15 minute ensemble time. Burst sampling was used to prevent aliasing of surface wave 
velocities. Sampling the first 40 seconds enabled the ADCP to measure the entire surface 
wave/swell cycle. By taking the average of this sample, the ADCP captured the true 


current velocity without significant bias caused by surface waves. 


Echo intensity is scaled, using about 0.45 decibels per WorkHorse count. It is a 
measure of signal strength. For the upward-looking ADCP, echo intensity 1s important. 
Echo intensity and the signal-to-noise ratio decrease with range. With increasing range 
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the echo intensity signal begins to weaken and flatten out. If a rapid increase (spike) in 
the echo intensity is observed, it is most likely due to the transducer’s signal side lobe 
having been reflected by the ocean surface and that the side lobes of the acoustic pulse 
are contaminating the velocity signal. Two criteria were used to flag bad data from the 
ADCP. For this study, once the echo intensity flattened out for several successive bins 
and then a spike in intensity occurred, it was then safe to assume that the ADCP stopped 
collecting good data. By comparing the echo intensity to the correlation magnitude, if 
both values fall below 64 counts then the data are flagged as bad. The second criterion is 
based on the three and four beam solution. If the percent good values of the 3 and 4 


beam solutions were less than 25% then the data are flagged as bad. 


In order to have one continuous year long time-series, SC 2 and SC 3 were 
patched together. There was a 2 hr 45 min. gap in time between deployments. The gap 
was first filled with NaN’s (not-a-number). For the spectral analysis, these were replaced 
with zeros. Time-series used for spectral analysis can not have NaNs or gaps. The bin 
depth values for each deployment were compared and matched with the closest single 
array of common mutual depth. The center of the lowest bin that was common to both 
deployments was 118 m and highest common usable bin was at 26 m. There were a total 


of 24 bins. 


The SBE 37-IM MicroCAT recorded conductivity, pressure and temperature. 
Temperature was measured by applying an alternating current to an airtight VISHAY 
reference resistor and an ultra-stable aged thermistor. Salinity was a derived from the 
conductivity and temperature measurements. Each mooring had three MicroCATs that 
sampled every 5 minutes. For the 1-year period the sampling depths were set to 50 m, 91 
m and 129 m. Portions of the conductivity data appeared to be offset, probably due too 
material in the cell; salinity for these periods were estimated using the T-S relationship as 


determined by good data. 


B. WINDS 


Wind data were downloaded from “NOAA’s National Data Buoy Center” 
(NDBC), http://www.ndbc.noaa.gov. Historical wind data were used from Stations 
46086, (San Clemente Basin, SCB), 46047, (Tanner Banks, TB), and 46025, (Santa 
Monica Basin, SMB). The period of time used was the same as the time-series. 
However, there was a data gap from May 10, 2007, until June 18, 2007, due to a failure in 
the anemometer at the SCB buoy (see Figure I-1 for reference location of the buoys in 


relation to the SCI). 
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Figure II-1. | Schematic diagram of the mooring at Eel Point. The same configuration 
was used for both SCORE 2 (SC2) and SCORE 3 (SC3). SC2 bottom depth: 130 m and 
SC3 bottom depth: 135 m. 
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Hl. DATA ANALYSIS METHODS 


A. INITIAL PROCESSING 


The ocean currents, temperature and salinity data collected from the ADCPs and 
MicroCATs for Eel Point were broken into different files for separate processing and 
clean up. Raw ADCP data were processed using the MATLAB® program 
Workhorseprol (written by Prof. Curt Collins) designed to flag and remove bad data 
points and interpolate with replacement values. The resulting time-series was then 
despiked and a linear interpolation was again used to replace the extreme values. These 
cleaned up interpolated data will henceforth be referred to as the raw data. Rotated, 
filtered, and vertically averaged data files were created for future use. The filter used to 
smooth the currents, and stratification (T-S) data was a 4" order low-pass Butterworth 
filter, the nondimensional cutoff frequency used for the currents was .01 and the cutoff 
for the T-S was .0035 (the cutoff values are different in order to compensate for the 
different sample rates), these values equate to an ~50 hour low-pass filter (SOlp). This 
was used to remove the tidal frequencies from the data. All data analyses were processed 


using tailored MATLAB® programs. 


B. UNSMOOTHED CURRENT OBSERVATIONS 


It was necessary to first analyze the raw data at all depths, in order to gain an idea 
of the vertical structure of the flow. Table 1 contains the mean, maximum, minimum and 
the standard deviation of the current velocity for the U and V components at each bin 
depth. The mean current flow was a dominant poleward flow in the V component and a 
dominant westward flow in the U component. The maximum (minimum) V component 
was recorded at 46.7 cm/s (-80.2 cm/s) and the U component was 61.1 cm/s (-58.4 cm/s). 

Histograms of speed and direction for the (unsmoothed) raw currents recorded by 
the ADCP at a 15 minute sample intervals for the near surface (34 m), middle water 


column (74 m), near bottom, and vertically averaged currents are shown in Figures III-1. 
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The speed and direction bin sizes used in the histogram were 0.5 cm/s and 1°. At the near 
surface, the most frequently observed current speed was 7.5 cm/s with a maximum 
current speed of 50 cm/s and most frequently observed direction of 333° T. In the 
middle water column, the most frequently observed current speed was 5.5 cm/s with a 
maximum current speed of 50 cm/s and most frequently observed direction of 337° T. In 
the near bottom currents, the most frequently observed current speed was 5.5 cm/s with a 
maximum speed of 44 cm/s and most frequently observed direction of 342° T. The most 
frequently observed current speed for the vertically averaged flow was 4 cm/s with a 
maximum speed of 37 cm/s and most frequently observed direction of 338° T. The 
histograms showed the predominant flow to be northward along the isobaths. The local 


isobaths were oriented at 330° T. 
C. PRINCIPAL AXES 


The ADCP data were rotated into their principal axes. The rotation maximizes 
the variance along the major axis and minimizes it along the minor axis. The major axis 
is now referred to as the alongshore component and the minor axis is the cross-shore 
component. The orientation of the principal axis, @ is determined by (Emery and 


Thomson, 1997): 


Where uw’ =u-—uandv'=v—v. Inthe new coordinate system, the rotated velocity 
components are defined as 
u.=u'cos@+v'sind 
and 
v. =—u'sind+v'cosé. 
In this coordinate system, the maximum variance is associated with v. which represents 


the alongshore (major axis) component, while the orthogonal component (minimum 


variance) is associated with u. the cross-shore (minor axis) component. The purpose of 
the principal axes were to determine the orientation of the variability of the flow and 
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decompose the velocity into uncorrelated scalar components (Emery and Thompson, 
1997). It also gives a statistical picture of the variance for a set of vector time-series 
(Table 2). 

For this study the principal axis rotation was calculated using a MATLAB® 
function, princax (written by Rich Signell). This function calculates @ (angle of 
maximum variance) in mathematical coordinates (90 is north and 0 east). An example 
for the vertically integrated flow is shown in a scatter plot (Figure I[II-2). For the 
example, @ = -55° (-55 subtracted from 90 gives 145 + 180 = 325). The alongshore 
direction 1s now orientated along 325° T to 145° T, and cross-shore 235° T to 55° T. The 
alongshore orientation aligned very closely with the topography. For example, NOAA 


chart 18762 indicated that the isobaths were oriented along 330 ° at the mooring location. 


D. SEASONAL VARIABILITY 


The bight region 1s a complex environment which has high temporal and spatial 
variability due to the local winds in addition to the eight nearshore islands, submarine 
canyons, ridges, and basins which disrupt the large scale circulation. This section 
describes of the seasonal variability of currents (vertically averaged over 24 bin depths 
between 26 m and 118 m) and temperature and salinity at the Eel Point mooring. The 
current velocity, temperature and salinity data used were divided into 3 month seasonal 
segments, starting with summer (July 10, — September 30, 2006), fall (October 1, - 
December 31, 2006), winter (January 1, - March 31, 2007) and spring (April 1, - June 30, 
2007). The data used to create all figures in this section were rotated and lowpass filtered 


(previously described) to remove the tidal variability. 
1. Summer 


In the summer, the vertically averaged currents (Figure III-3) show fairly weak 
alongshore currents which are predominantly poleward flowing. However there are 3 
short current reversals. The cross-shore currents are very weak and flow toward the island 
until the end of August when a significant increase in current velocity (15 cm/s) along 
with a flow reversal offshore occurs for about 15 days. The contour plots (Figures III-4) 
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show the current velocities with depth. The contours are consistent with the vertically 
averaged currents and show that the cross-shore reversals penetrate throughout the water 


column, but were stronger at the surface and weaken toward the bottom. 


The summer temperature time-series (Figure III-5) shows a significant difference 
in temperature between 50 m and 91 m, AT averaging 3° C between these depths. At 50 
m, the temperature steadily rose to a maximum of 14.3° C and maximum AT of 4.5° C 
between depths 50 m and 129 m in mid-September. Temperatures at 91 m and 129 m 
increased little. For summer salinity (Figure HI-5), the salinity increased with depth 
compared to temperature which decreased with depth. This figure shows that salinity and 


temperature were inversely proportional to each other. 
2 Fall 


The vertically averaged currents in the fall (Figure III-6) start off weak with 3 
small reversals in current flow in the alongshore and cross-shore flow. By mid- 
November the alongshore flow increases and is predominantly poleward in direction and 
the cross-shore flow increases eastward (toward the island). In the fall, contour plots 
(Figure III-7) of the alongshore flow and cross-shore flow almost mirror each other. 


When the alongshore flow is poleward the cross-shore flow is eastward. 


The fall temperatures (Figure III-8) were slightly warmer than the summer at 50 
m with a maximum peak 15.1° C and maximum AT 5° C at end of October. After the 
maximum peak, the temperatures decreased and averaged between 13-14° C. There was 
a significant decrease in AT= 2.2” C between 50 m and 91 m. There were variations at 50 
m which did not propagate down through the water column. There were no significant 
variations 1n the salinity time-series (Figure III-8). The relative stratification between the 
layers, remain constant, until late December when more mixing was observed between 


the layers. 
3. Winter 


The winter vertically averaged currents (Figure III-9) show a fairly strong and 


predominantly alongshore poleward flow. At the beginning of February, there was an 
16 


onset of a strong equatorward reversal lasting for 10 days. The cross-shore flow 
fluctuated at the beginning of January. It then remained relatively weak throughout the 
season. The winter alongshore contour plot (Figure [II-10) shows that the oscillations 
and reversal observed in the vertically averaged currents do penetrate down through the 
water column, with maximum velocities observed at mid-depths. In the cross-shore the 


eastward flow was uniform throughout the water column. 


The winter temperatures (Figure III-11) were the warmest at all depths in 
comparison with the other seasons. At 50 m there were a lot of variations in temperature. 
The maximum temperature peak was 14.8° C with a max AT 4.8° C between the depths. 
At 91 m and 129 m the temperatures were 1-2° C warmer than summer and fall. The 
winter salinity time-series (Figure IHI-11) shows less stratification between the 50 m and 


91 m, meaning there was more mixing between the two depths. 
4. Spring 


The spring vertically averaged currents (Figure III-12) show the strongest (23 
cm/s) poleward alongshore flow. The cross-shore flow was relatively weak and directed 
eastward. Starting around April 10™ the alongshore flow showed an oscillation with a 2-3 
day period which lasted about 10 days (hereafter, referred to as the April event). The 
oscillation can be seen in the contour plot (Figure [JJ-13): the current velocities were 
strongest at mid-depth. The vertical phase change for the April event was estimated to 


be 0.12 cm/s, (Az = 40 m determined for a 2 day period (April 15 — 17, 2007). 


In the spring temperature and salinity time-series (Figure I-14), the April event 
(oscillation) was observed at all depths. The temperature steadily decreased with time at 
all depths. The largest temperature difference was observed at 50 m, decreasing from 
14.1° C to 11.4° C (AT=2.7° C). At 91 m and 129 m the temperature steadily decreased 
to < 10° C. The salinity stayed relatively stratified throughout the water column until 
early June when it increased, weakening the stratification. The observed cooler 


temperatures and saltier water characteristics are mostly likely the result of upwelling. 
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>; Recap of Seasonal Variation for the Year 


Looking at the entire time-series, there were many 2-3 and 3-4 day patterns of 
oscillations observed in the current velocity and stratification. However, there were no 
discernable patterns over the whole year. Summer had the weakest observed currents. 
Winter and spring had the strongest currents. Winter had the warmest temperatures and 
spring had a cooling trend along with the coldest observed temperatures. Stratification 


was strongest throughout the summer and weakest at the end of spring due to upwelling. 


E. SPECTRAL ANALYSIS 


Spectral analysis is a way to show how the variance of a time-series varies as a 
function of frequency. Three types of spectra were calculated for this thesis, the kinetic 
energy, rotary, and wavelet spectra. The kinetic, rotary and wavelet spectra were 
computed using MATLAB® programs. The MATLAB® programs were based on the 
discrete fast Fourier transform (FFT) and the periodogram method. A Hanning window 
was applied with a 50% overlap. Prior to computing the spectrum, the record mean and 
trend were removed from the time-series. A 180 day piece length and a 95% confidence 


interval were used for both kinetic and rotary spectra. 
1. Kinetic Energy Spectra 


Kinetic energy (KE) spectra are a useful tool to study deterministic and stochastic 
processes. The Eel Point current time-series contained a mixture of both deterministic 
and stochastic processes. Deterministic processes are predictable, such as the tides, while 
stochastic processes are random, such as internal/inertial waves. The raw data were used 
to construct the KE spectra. The KE spectra were used specifically in this thesis to show 
the frequency of the dominant energy. Variance-preserving KE spectra were generated 
for three levels within the water column; 110 m, 74 m, and 34 m (Figure III-15). It 
showed that semi-diurnal and diurnal tides dominated the energy spectra. The semi- 
diurnal tidal energy was strongest at mid-water and weakest at the near-bottom. The 


diurnal tide was strongest at the near-surface and near-bottom. Knowing that the 


18 


dominant energy was in the semi-diurnal and diurnal frequencies, a band-averaged KE 
spectra at the same depths was created (Figure III-16), it confirms that the diurnal and 
semi-diurnal KE dominated the spectrum and that the semi-diurnal KE was about 20% 


greater than the diurnal KE. 
2. Rotary Spectra 


Rotary Spectra are used to separate the velocity vector frequency, @, into 
clockwise and counterclockwise rotating circular components. The circular components 
replace the Cartesian components (u, v). The vector addition of the two oppositely 
rotating circular components cause the tip of the combined vector to trace out an ellipse 
over one complete cycle. If one of the two components equals zero, the motion is 
circular polarized, and if the two components are equal (same magnitude), then the 
motion 1s rectilinear (back and forth). In most cases, rotary spectral analysis has proven 
to be useful when investigating current variability. For the purpose of this thesis, rotary 
spectra were used to investigate inertial motions, diurnal frequency and any other types of 
oscillatory flow at the near-bottom (110 m), mid-water column (74 m), and near surface 


currents (34 m). 


Figure III-17 shows that the clockwise energy was dominant over the 


counterclockwise energy at all three depths. The inertial frequency @=2Qsin@ was 


calculated using a latitude of 32°, resulting in the inertial period 1/@= ~22.5 hours. At 
the near surface and near bottom the dominant peaks were at the diurnal, inertial and 


semi-diurnal. The mid-depth, the spectra did not contain the inertial peak. 


F. TIDAL ANALYSIS OF CURRENTS 


Tidal variability is often the largest signal of energy in a typical oceanic time- 
series (Pawlowicz et al, 2002). Tides are deterministic, meaning that they can be 
predicted and are mostly dominated by diurnal and semi-diurnal frequencies. Coastal 
environments such as the Southern California Bight region are highly dynamic in nature. 
The currents are often affected by deep ocean variability such as the tides which 


propagate through the shallower waters as a wave or combination of waves. In order to 
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understand the nature of the tidal currents, a harmonic (tidal) analysis was used to obtain 
the tidal constituents. The tidal analysis allows specific frequencies to be examined and 
amplitudes and phases to be estimated in a least-square sense (Emery and Thomson, 


1997). 


For this thesis, the Forman tidal harmonic analysis program was run using the 
MATLAB® code from Pawlowicz et al., (2002). The analysis uses the MATLAB® 
program T TIDE to run the tidal analysis using the raw data. T TIDE uses the least- 
Squares method which is constrained to fit the frequencies that Forman specified as 
tidally important. The advantage of using the least-squares method is that it permits 
resolution of 67 tidal constituents (for this data set) of which 45 are typically 
astronomical in origin. Tables III-3 and III-4 list the five dominant constituents found 
using the tidal analysis. Further dynamical analysis requires the separation of the tidal 
effects from inertial effects. The semi-diurnal (M2) and diurnal (K1) tides were the 
dominant frequencies throughout the entire time-series, they were also very prominent in 
the KE spectra. The T TIDE program creates a data file which contains all tidal 
constituents found at Eel Point mooring. The vertically averaged tidal constituents were 
computed and used to compute the predicted tides using the MATLAB® program 
T PRED. The predicted tides were defined as the barotropic tidal currents. These 


barotropic tidal currents were used later in the wavelet analysis. 


G. WAVELET ANALYSIS 


Wavelet analysis is a good tool for localizing time-frequency information without 
requiring that the time-series be stationary. Wavelet analysis involves transforming a 
time-series into a series of coefficients that contain all the information of the original 
signal, which is similar to Fourier analysis. In a wavelet analysis, the signal is 
reconstructed by adding wavelets of different scales and at different translations together; 
as in a similar manner Fourier analysis reconstructs by adding sines and cosines of 
different periods together. The wavelet analysis optimally tiles the time-frequency 
domain allowing for dominant modes of variability to be determined and to distinguish 


how those modes vary in time. 
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For this analysis, the methods described in Torrence and Campo (1998) “Guide to 
Wavelet Analysis” and their toolbox were used to generate the wavelet spectra. Prior to 
computing the wavelet analysis, the record mean, trend, and barotropic tidal currents 
were removed from the time-series. A continuous wavelet transform which can use both 
orthogonal or nonorthogonal wavelet functions turns a signal into a function with two 
variables, scale and time. The basis function used was a Morlet wavelet which is the 


product of a complex exponential plane wave modulated by a Gaussian: 


-1/4_ ian 7/2 
9 


Wi(7)=a ee 


where 7 is the nondimensional “time” parameter and @, is the nondimensional 


frequency. In a wavelet power spectrum, different modes of energy are contoured in 
order to interpret how the amplitude of the signal changed with time and frequency. The 
wavelet function is normalized at each scale s to have unit energy. This unit energy 
criteria ensures that the wavelet transforms at each scale s are directly comparable to each 


other and other time-series. The normalization process takes place by using: 


[ee -(2) {ee 
S S S 


where s is the wavelet scale and n is the translation parameter. The factor of s "7 


indicates y,(7) has been normalized. 


The Morlet wavelet was used for this analysis since it is particularly good at 
isolating periodic signals. The width of the wavelet function is defined in Torrence and 
Compo (1998) as the e-folding time of the wavelet amplitude. The resolution 1s 
determined by the wavelet scales. Small wavelet scales resolve the high frequency 
components, the wavelet itself is broad in frequency and the amplitudes of the peaks in 
the spectrum are smoothed out. The wavelets are scaled based on a relationship between 
the equivalent Fourier period. Because wavelets deal with finite length time-series, error 
will occur at the beginning and end of the wavelet power spectrum. One solution is to 
pad the end of the time-series with zeroes before computing the wavelet transform and 


then remove them afterward (Torrence and Compo, 1988). The cone of influence (COI) 
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is the area of the wavelet spectrum where edge effects become important. Outside the 


COI, estimates contain errors and may not represent the true period or energy. 


For this analysis, Torrence and Compo (1998) software was used with the time- 
series padded with zeros. Ten suboctaves per octave (an octave is the difference between 
2 frequencies that relate to each other by a factor of 2) were used to characterize the 
resolution in frequency. The highest resolved frequency corresponded to a period of 30 
minutes and the lowest resolved frequency to 44 days for SC2 and SC3. The lowest 


resolved frequency for the 2 month segment was 20 days. 


The barotropic tidal currents were subtracted from the rotated raw data to create 
the baroclinic tidal currents. The baroclinic tidal currents and the wind stress (discussed 


later) thus created were used in this wavelet analysis. 


The wavelet analysis will look at SC2 and SC3 separately and then look at a 
smaller two month time segment. The same three depths used in the previous spectral 
analysis (34 m, 74, 110m) will also be the main focus of this analysis. The wavelet 
spectra depth figures represent the near-inertial/diurnal energy. They were created by 
averaging three different periods: 21.8, 23.4, and 25.1 hours which are not necessarily 
independent, estimates of wavelet power. SC2 alongshore and cross-shore wavelet 
components with depth (Figure III-18) show that the energy was only at the near surface 
and near bottom depths. In the middle of the water column (60-70 m) there was very 
little energy. SC2 (Figures III-19, 20, & 21) show that at all three depths the there was 
higher energy at the diurnal and semi-diurnal periods, the energy was aperiodic. A higher 
energy band at the 6 — 10 day period was present at all three depths (circled in red), being 


strongest at 74 m (this will be compared later with to the wind stress wavelets). 


In the SC3 alongshore and cross-shore wavelet spectra with depth (Figure III-22) 
the same results were obtained as in SC2. The energy occurred in the near surface and 
near bottom with very little in the middle of the water column. At all three depths in SC3 


(Figures HI-23, 24, & 25) there was more dominant diurnal energy in the alongshore flow 
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and semi-diurnal current in the cross-shore flow. Between April 12 — 21, the April event, 
there was strong 2-3 day period with energy present all three depths in the alongshore 
flow and strongest at 74 m (Figure II-24). 


In order to investigate the April event in greater detail, a two month (March 15 — 
May 15, 2007) period was analyzed. The rotary wavelet power spectra for all depths 
(Figure III-26) were computed in order to take a closer look at the diurnal-inertial band. 
Clockwise energy was dominant over the counterclockwise energy. In the clockwise 
power spectrum there was higher energy in the upper layer and bottom layer, however the 
energy peaks between the two layers were not contiguous. At ~65 m there was little 
energy present. Looking at the rotary wavelet spectra at the same three depths as before, 
(Figure III-27, 28, & 29) there were intermittent periods of enhanced near-inertial/diurnal 
and semi-diurnal energy in the clockwise and very little seen in the counterclockwise 
spectra. Between April 1 and April 6 there was a very strong event at 34 m that was 
dominated by the clockwise energy but still fairly strong in the counterclockwise (circled 
in red in Figure HI-27). Nearly rectilinear motion only appeared at the near surface (34 
m) depth. The April event was not present near surface. The clockwise wavelet spectra 
had high energy at 2 days decreasing to | day (near inertial/diurnal period) and the 
counterclockwise wavelet spectra had very weak energy during the event time frame. At 
mid-depth (74 m), the clockwise and counterclockwise energy appeared equal during the 
April, making it rectilinear (Figure HI-28). The near bottom (110 m) counterclockwise 
energy was dominant over the clockwise energy, meaning the currents were slowly 
rotating in the counterclockwise direction (Figure I[I-29). At all three depths 
intermittent modes of higher energy at the near-inertial/diurnal and semi-diurnal periods 
(only in the clockwise rotary wavelet) were observed. The strongest intermittent modes 
of near inertial/diurnal and semi-diurnal energy were observed at the near bottom depth. 
Over all, the clockwise energy was dominant at all depths. The near surface alongshore 
and cross-shore velocity components (Figure III-30) show the same high energy mode at 
the near-inertial/diurnal period that was coincident with the rotary wavelet spectra. The 


April event was very weak in the near surface in the alongshore flow. However, at the 
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mid-depth (Figure HI-31) the April event was strongest in the alongshore component, and 


it also appears at near bottom (Figure III-32). 


Wavelet power spectra were computed from the U and V wind stress components 
for all three buoy stations: TB, SMB, and SCB. The time frame for the wind stress 
spectra ran from July 10, 2006 through May 10, 2007 (the data gap at SCB buoy, 
prompted this ending time). Little diurnal and semidiurnal energy was observed at the 
three wind stations. Most of the energy was observed at periods greater than two days. 
Higher energy observed at the 6 — 10 day period (Figure’s II- 33, 34, & 35) beginning 
September 6, 2007 through the beginning of October 2007, was coincident with the 
higher energy observed in the currents for the same time period. This higher energy 
could be due to a storm passing through the area. Based upon these observations, the 
April event does not appear to be induced by the wind stress from any of the three buoy 


locations. 


H. WIND STRESS FORCING 


1. Wind Stress and Curl 


Many large and small scale ocean circulations are dependent upon the interaction 
of the atmosphere and ocean. Winds in the bight were generally weak, but were highly 
variable in comparison to winds north of Point Conception. The wind stress values (the 
same values which were used in the wavelet analysis) for all three NDBC buoy stations: 
TB, SCB, and SMB were computed using the procedure suggested by Large and Pond 
(1981) using a neutral drag coefficient. The wind stress time-series for all three stations 
(Figure III-36) were computed to compare the sheltering effects at each station. SCB and 
SMB have relatively weak wind stress values. Tanner Banks is on the western edge of the 
bight and the open ocean. The TB buoy was not affected by the sheltering that occurred 


within the bight so naturally the wind stress was going to be much larger at this location. 


The curl of the wind stress was calculated from these three NDBC buoy stations. 
The wind stress and the alongshore currents were smoothed to 10 days using a 4°" Order 


Butterworth filter. The curl of the wind stress was compared to the alongshore currents 
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(Figure III-37). This was done to investigate the possible forcing of the current reversals. 
No visual correlation was found and a cross-correlation computation confirmed this 


conclusion (figure not shown). 
2. Wind Forced Slab Model 


The location of Eel Point mooring, 32° N latitude, is considered a critical latitude 
for diurnal-inertial resonance. At 30° latitude, the inertial frequency is equal to the 
frequency of the diurnal tide (Craig, 1989). In order to investigate if inertial oscillations 
are prevalent at this location, a simple slab model following Pollard and Millard (1970) 
was created. The simple slab model uses the wind stress to predict if there are any 


inertial oscillations as result of the forcing. 


The wind stress calculated from the SCB buoy was used to force the model of 
predicted inertial currents. The model uses a simple slab for input of the stress from the 
atmosphere. This model follows that discussed in Paduan et al. (1989), using the 


formulation from Pollard and Millard (1970). In order to compute the inertial currents, 
Q=u+iv (Q 1s the inertial current) is driven by the surface wind stress pY =7* +ir°. 


The wind stress is taken up uniformly over a slab of depth h (which is treated as a 


constant): 


dO Y 
EO Se. 
aa 


(equation | in Paduan et al. 1989) where (2 =r + if and r is the Raleigh frictional decay 
included with a an e-folding time of r'= 88 hours. The solution is: 


iy. 
_ —QOr —Q(t-t') 
O=Q,e" + |, Py dt , 


(equation 2 in Paduan et al. 1989) where Q, is the initial current and equal to zero for this 


study. 


This model uses complex demodulation at the near inertial frequency. It removes low 


frequency and tidal motions which are not expected to be correlated with the local wind. 
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For this analysis, unsmoothed wind stress values computed from SCB buoy were 
used in order for the model to predict oscillatory current responses. The mixing layer 
was set to h = 34 m for comparison with the complex demodulated near surface currents 
at depths 34 m. There were periods of inertial motion output with this model. The 
currents at 34 m do not reflect any inertial frequency oscillations predicted by the model 
(Figures not shown). However, there was a strong near inertial/diurnal mode of energy 
in the 2-month wavelet spectra in the along and cross-shore flow. This higher energy was 
only observed at the near surface. The demodulated currents at 34 m showed a peak in 
the current speed at the same time of the observed higher energy in the spectra (around 
April 4, 2007), the inertial oscillation was not predicted by the model using the SCB buoy 
wind stress. The model was run using the demodulated currents at 38 m and there was no 
peak in the current speed observed, so the energy was only associated with near surface 
currents < 38 m. It was possible there may have been other inertial frequency oscillations 
in the near surface which did not penetrate to 34 m. These could not be determined with 
this data set, since the top bin was so deep and wind stress used in the model was not 


localized to the Eel Point mooring. 
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Table 1. 


ADCP statistical results from the observed (unsmoothed/unrotated) data. U and V 


component values of the mean, maximum, minimum and standard deviation. 


Standard 


Deviation 


U (cm/s) 





Standard 
Deviation 


V (cm/s) 


Depth Mean U Max U Min U Mean V Max V Min V 
(cm/s) (cm/s) (cm/s) (cm/s) (cm/s (cm/s) 


-28.9 2.8 38.6 
ee! 


a 
13 
i 


34m -1.8 27.9 -31.3 2.9 42.0 -21.0 

38m 28.0 -30.2 41.8 299 

42 m 28.5 -29.8 44.0 -22.8 6.2 

46m 28.4 -32.3 44.2 -24.9 (Pa) 
50m 29.4 -33.2 40.6 -29.1 
58m 28.6 -33.0 39.0 -27.3 | “6. 

102 m 293 28.9 -42.4 4.4 44.1 -27.0 89 
106 m 32.2 AD 9 44.8 -32.9 8.2 
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Table 2. —- Principal Axes. 


Semi-Major Axis 


Direction Semi-Major Axis Semi-Minor Axis 
Magnitude (cm/s) Magnitude (cm/s) 
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Table 3. |The 5 dominant tidal constituents computed from the vertically averaged currents. 
Ol and K1 are diurnal and N2, M2 and S2 are semi-diurnal tides. 


Period Major Axis Minor Axis 
Tidal Constituents eal a Sl Greenwich Phase (deg) 


| 580 82 
23.93 
12.66 
12.42 


=. 7 et ge ae ae 





Table 4. The 5 dominant tidal constituents computed from the 50 m MicroCAT pressure 
data. 


Tidal Constituents Amplitude (cm) Greenwich Phase (deg) 
218.65 
C K™~dSSC“‘“CSWC 207.14 


a > 163.30 
(70.41 
142.10 
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Figure III-1. Histogram of speed (left) and direction (right) for ADCP currents, at near 
surface (34 m), mid-water (74 m), near bottom (110 m), and vertically averaged. 
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Figure III-2 Scatterplot generated from the Principal Axis Rotation. @ = 325° T. 
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Figure III-3. Summer vertically averaged currents (data were rotated and vertically 
averaged over 24 bin depths between 26 m— 118 m). 
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Figure III-4. Summer current velocities for the alongshore and cross—shore flow with 
depth. Alongshore flow was contoured from 35 to -35 at 5 cm/s intervals and the cross- 
shore flow was contoured from 15 to -15 at 2 cm/s intervals. 
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Figure III-5. Summer temperature and salinity measurements collected from the 
MicroCATs at depths: 50 m (red), 91m (blue), and 129 m (green). 
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Figure III-6. Fall vertically averaged currents (data were rotated and vertically averaged 
over 24 bin depths between 26 m— 118 m). 
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Figure III-7. Fall current velocities for the alongshore and cross—shore flow with depth. 
Alongshore flow was contoured from 35 to -35 at 5 cm/s intervals and the cross-shore 
flow was contoured from 15 to -15 at 2 cm/s intervals. 
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Figure III-8. Fall temperature and salinity measurements collected from the 
MicroCATs at depths: 50 m (red), 91m (blue), and 129 m (green). 
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Figure III-9. Winter vertically averaged currents (data were rotated and vertically 
averaged over 24 bin depths between 26 m— 118 m). 
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Figure IJJ-10. Winter current velocities for the alongshore and cross—shore flow with 
depth. Alongshore flow was contoured from 35 to -35 at 5 cm/s intervals and the cross- 
shore flow was contoured from 15 to -15 at 2 cm/s intervals. 


39 


—— 40m 
——9/m 
aa {9 m 


VYWWAL 


e 
z 
E 
2 


late 251 tro) sen Dl rl 


ail 
Janay | = Mareh 31, TNF 





Figure III-11. Winter temperature and salinity measurements collected from the 
MicroCATs at depths: 50 m (red), 91m (blue), and 129 m (green). 
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Figure III-12. Spring vertically averaged currents (data were rotated and vertically 
averaged over 24 bin depths between 26 m— 118 m). 


4] 


ay 4F ld 42) 6438 ee sh oT) | 6ae 


S? Sid Sz) 828 el 
April 1} - June 30, 2007 


(ross-shore Flow 


a4 


2) S28 Gl wf S14 GI 
April | - June 30, 2007 





Figure IJIJ-13. Spring current velocities for the alongshore and cross—shore flow with 
depth. The oscillation is evident from the 12™to the 21“ of April, which is referred to as 
the April event and circled in red above. (contour intervals are as previously mentioned). 
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Figure III-14. Spring temperature and salinity measurements collected from the 
MicroCATs at depths: 50 m (red), 91m (blue), and 129 m (green). 
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Kinetic Energy Spectra at Depth 34 m 
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Figure IJI-15. Variance Preserving Kinetic Energy Spectra for near-surface (34 m), 
middle water column (74 m), and near-bottom (110m). The spectra were dominated 
by energy in the semi-diurnal (M2 peak) and the diurnal (K1) tidal frequencies. 
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Figure IIJ-16. Band averaged KE spectra, showing that dominant energy 1s 1n the semi- 
diurnal and diurnal frequencies. 
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Figure IJJ-17. Rotary Spectra for near-surface (34 m), mid-water (74 m) and near bottom 
(110 m). Red (green) is clockwise (counterclockwise). 
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Figure III-18. SC2 wavelet spectra of the near-inertial/diurnal band with depth for the 
alongshore (v) and cross-shore (u) components. (Contour intervals are non-linear) 


47 


Wavelet Power Spectrum for the Alongshore velocity ( 


Period (days) 


01,0625 4 


7/30 08/19 Oo8 9/28 LOLs 1 Laer L127 [2/17 0 Lk 


Wavelet Power Spectrum for the Cross-shore velocity (u) at z=34 m 


Period (days) 


OF30 Og/L9 STAD. 09/28 LOLs L LA? 





Figure III- 19 SC2 wavelet power spectra from July 10, 2006 — January 18, 2007, depth 
34m. Dashed white line is the COI. Higher modes of energy at the 6-10 day period, 
circled in red. (Contour intervals are non-linear) 
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Wavelet Power Spectrum for the Alongshore velocity (v) at z=74 m 
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Figure IJI-20. SC2 wavelet power spectra from July 10, 2006 — January 18, 2007, depth 
74 m. Higher modes of energy at the 6-10 day period are circled in red. 
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Figure III-21. SC2 wavelet power spectra from July 10, 2006 — January 18, 2007, depth 
110 m. Higher modes of energy at the 6-10 day period are circled in red. 
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Figure III-22. SC3 wavelet spectra of the near-inertial/diurnal band with depth for the 
alongshore (v) and cross-shore (u) components. 
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Figure IJ-23. SC3 wavelet power spectra from January 18, 2007, - July 23, 2007 depth 
34 m. Dashed white line is the COI. 
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Figure III-24. SC3 wavelet power spectra from January 18, 2007, - July 23, 2007 depth 
74m. Higher mode of energy at the 2-3 day period are circled in red (April event). 
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Wavelet Power Spectrum for the Alongshore velocity (v) at z=110 m 
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Figure IIJI-25. SC3 wavelet power spectra from January 18, 2007, - July 23, 2007 depth 
110 m. Higher mode of energy at the 2-3 day period are circled in red (April event). 
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Figure IIJ-26. Rotary wavelet power spectra for the period: March 15, - May 15, 2007, 
clockwise (counterclockwise) is on the right (left). 
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Figure IJJ-27. Rotary wavelet power spectra for the period: March 15, - May 15, 2007, 
depth 34 m. Higher modes of energy at the near inertial/diurnal period and at the 2-1 day 
period are circled in red. 
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Figure IIJ-28. Rotary wavelet power spectra for the period: March 15, - May 15, 2007, 
depth 74 m. Higher modes of energy at the 2-3 day period are circled in red (April 
event). 
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Figure IJJ-29. Rotary wavelet power spectra for the period: March 15, - May 15, 2007, 
depth 110 m. Higher mode of energy at the 2-3 day period are circled in red (April 
event). 
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Wavelet Power Spectrum for the Alongshore velocity (v) at z=34 m 
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Figure IJI-30. Wavelet power spectra for the period: March 15 — May 15, 2007 at 34 m. 
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Wavelet Power Spectrum for the Alongshore velocity (v) at z=74 m 
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Figure IJJ-31. Wavelet power spectra for the period: March 15 — May 15, 2007 at 74 m. 
Higher mode of energy at the 2-3 day period are circled in red (April event). 
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Figure IIJ-32. Wavelet power spectra for the period: March 15 — May 15, 2007 at 110 
m. Higher mode of energy at the 2-3 day period are circled in red (April event). 
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Figure III-33. Wind Stress wavelet power spectra for winds from NDBC - San Clemente 
Basin (SCB) buoy. Higher mode of energy at the 6-10 day period are circled in red. 


62 


Wavelet power spectrum for U (SCB) 


- bee Re a 


ONT GHZ MEK ARG 


PEG CR OIG TAR GEE WIE ILS? LF MIT GL! OL Is 





Figure III-34. Wind Stress wavelet power spectrum for winds from NDBC - Santa 
Monica Basin (SMB) buoy. Higher mode of energy at the 6-10 day period are circled in 
red. 
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Figure III-35. Wind Stress wavelet power spectra for winds from NDBC — Tanner Banks 
(TB) buoy. Higher mode of energy at the 6-10 day period are circled in red. 
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Figure IJJ-36. Wind stress time-series for NDBC stations: 46025 — Santa Monica Basin, 
46086 — San Clement Basin, and 46047 — Tanner Banks. 
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Figure IJJ-37. Wind stress curl calculated between TB, SCB, and SMB NDBC stations, 
compared to the smoothed vertically averaged alongshore currents. 
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IV. SUMMARY OF RESULTS 


This study was the first to describe the year-round current structure and 
stratification at Eel Point and is part of a continuing study. The currents around SCI 
were energetic and mostly dominated by semi-diurnal and diurnal tides. Periods of 
higher energy were observed, perhaps the most interesting of which was a ten day period 
in April, 2007. This chapter will summarize the important results from the study and 
place them in the context of other regional studies. Ocean currents and tides, seasonal 


variability of stratification, and forcing by wind stress and curl will be discussed. 


A. OCEAN CURRENTS AT EEL POINT MOORING 


Both the unsmoothed and smoothed currents were dominated by poleward flow 
along the isobaths with a mean (max) speed of 4 (50) cm/s. Occasional current reversals 
occurred that did not appear to be seasonal; on average 2-3 reversals occurred per season 
with the strongest reversal in the winter with observed equatorward maximum speed of 
30 cm/s. In the summer the currents were relatively weak, with one significant (10 day) 
reversal in the cross-shore flow. In fall, currents remained relatively weak until mid- 
November when the poleward alongshore flow increased as did the eastward cross-shore 
flow. The winter currents were highly variable; alongshore currents were more strongly 
poleward than in the fall with the largest flow reversal at the beginning of February. The 
spring season had the strongest alongshore currents while the cross-shore flow was weak, 


with a mean flow of 1-2 cm/s. 


There was only one other record of direct current measurements in San Nicholas 
basin for comparison with these results. In 1989, Tsai et al., (1990) measured the 
submesoscale current structure in San Nicholas Basin using synoptic shipboard current 
surveys during July 17-21 and September 2-6; the strongest currents they observed were 
poleward alongshore currents (mean flow 15 cm/s) on the eastern side of the basin, 
centered at 40-60 km to the west of SCI. These observations support the direction and 
magnitude of the currents reported here at Eel Point. 
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Next, seasonal variability of currents off Eel Point was compared to seasonal 
current variability observed at about the same latitude on the continental shelf offshore of 
Del Mar, California. The Del Mar observations were 3.6 km from shore in water about 
65 m deep versus 5 km offshore of SCI in 135 m depth at Eel Point. The Del Mar current 
observations were made in the summer (July 27 to September 11, 1978), fall (October 21 
to December 4, 1978), and winter December 21, 1978 to March 27, 1979) (Lentz and 
Winant, 1986). For this comparison, the rotated unsmoothed data were used, sorted into 
the same three seasons as in the Del Mar study, and a depth of 60 m was chosen. 

Results are tabulated in Table 5. As expected, mean flow at both locations was 
along the isobaths. Weak poleward flow was observed at Del Mar during the summer 
(1.3 cm/s), weak equatorward flow observed in the fall (-1.1 cm/s), and a stronger 
equatorward flow observed in the winter (-6.8 cm/s). At Eel Point, the mean poleward 
flow was lowest in the summer (1.3 cm/s) and stronger in the fall (5.0 cm/s) and winter 
(5.9 cm/s). So, at 60 m, the currents were about the same at both locations in summer but 
as the year progressed, the current accelerated in opposite directions. 

The standard deviation of the 60 m flow at Del Mar and Eel Point were about the 
same 1n summer, ~ 6 cm/s, but increased in fall and winter to ~ 10 cm/s at Eel Point, but 
remained ~6 cm/s at Del Mar. The standard deviation of the across shelf flow was much 


smaller at Del Mar because of the proximity of the current observation to the bottom. 
B. SEASONAL VARIABILITY OF STRATIFICATION 


As is typical for Eastern Boundary regions, stratification was minimum in 
summer and maximum in October. In July, the minimum temperature difference between 
50 m and 129 m was observed (1.2° C) and maximum temperature difference between 
these depths was 4.5° C in October. The temperature change between 50 m and 91 m 
was about twice that observed between 91 m and 129 m. Salinity increased with depth so 
it added to the stratification. The maximum salinity difference between 50 m and 129 m 
was AS=0.5 and occurred in July 2006 and the minimum salinity difference was AS=0.3 
in April 2007. The salinity change was about the same between 50 and 91 m as between 


91 mand 129 m. 
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At all depths, maximum salinity and minimum temperatures were observed in 
July 2007, and minimum salinity in October 2006 and maximum temperature in February 
2007. The July events were probably related to seasonal upwelling and it was noted that 
the subsurface temperature minimums occurred about 3-4 months later at Eel Point than 
the surface temperature minimums at inshore locations, e.g. NOAA Buoy 46025 
(Catalina Ridge, CA) (Goericke et al., 2007). Minimum subsurface salinity in the region 
should be correlated with the onshore flow in the Southern California eddy but the 
observed October minimum was at odds with the fall to winter occurrence of the 
Southern California eddy. This suggests that pycnocline displacements, (as opposed to 
changes of water mass properties on isopycnals) control the temporal variability of 


salinity in San Nicholas basin. 
C. TIDAL RESULTS 


The dominant tidal energy was found in the semi-diurnal and diurnal tidal energy 
bands. The KE spectra and band averaged KE spectra showed that the dominant energy 
was associated with the semi-diurnal and diurnal tides at all three depths (34 m, 74 m, 
and 110 m). Semi-diurnal was strongest at mid-depth and diurnal was strongest at the 
near bottom and near surface. Semi-diurnal tidal peaks were dominant at all depths in the 
rotary spectra, followed by diurnal peaks. Inertial energy was observed at the near 


surface and near bottom depths in the rotary spectra. 


The unsmoothed vertically averaged currents were dominated by semi-diurnal and 
diurnal tides (barotropic tidal currents). Semi-diurnal and diurnal modes were the 
dominant energy observed in the baroclinic tidal currents. In the (2 month) wavelet 
power spectra (which used the baroclinic tidal currents), there were strong pulses of 
energy in the near inertial/diurnal along and cross-shore currents that was only observed 
at the near surface depth. Intermittent pulses of higher energy at the near-inertial/diurnal 
and semi-diurnal periods energy in the along and cross-shore and rotary wavelet spectra 
were also observed. The amplitude of these pulses increased with depth so the near 
bottom had the highest energy associated with the near inertial/diurnal and semi-diurnal 


currents. 
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The amplitude of the semi-diurnal baroclinic currents at Eel Point did not track 
with the spring/neap cycles. Lerczak, (2000) studied the internal waves on the 
continental shelf near Mission Beach, California located at 32.75° N latitude (Eel Point 
mooring is located at 32.92° N latitude), and also found that semi-diurnal currents on the 
shelf were intermittent and did not follow the spring/neap cycle of the barotropic tides. 
Lercezak, (2000) also found intermittent diurnal currents which were driven by a 
persistent sea breeze that occurred off the coast of southern California. The intermittent 
diurnal currents observed at Eel Point could not be linked to the winds observed at 


regional NOAA buoys. 
D. APRIL (OSCILLATION) EVENT 


In April 2007, the alongshore flow oscillated with a 2-3 day period for 
approximately 10 days. These oscillations propagated downward through the water 
column with a vertical phase velocity of ~0.12 cm/s and maximum alongshore velocity 
(30 cm/s) occurred at the mid-depths. This event was not observed in the cross-shore 
velocity. This oscillation was not observed near the surface (34 m). The clockwise 
wavelet spectrum had high energy at 2 days decreasing to 1 day (near inertial/diurnal 
period); at this time the counterclockwise wavelet spectrum had a very weak pulse of 
energy. At mid-depth (74 m), the highest energy was observed and the clockwise and 
counterclockwise energy were equal, confirming that the flow was rectilinear alongshore. 
Near the bottom (110 m), the counterclockwise energy was dominant over the clockwise 


energy, so the horizontal currents were rotating counterclockwise. 


During the period of these current oscillations in April, there was no energy 
observed in the wind stress wavelet spectra at any of the three buoy stations. The April 


event was therefore not induced by the regional wind stress. 


E. WIND STRESS AND CURL 


Wind stress within the bight 1s weak in comparison of wind stress north of Point 
Conception (Hickey, 1998). A 10-year time-series (Dorman and Winant, 1995; Winant 


and Dorman, 1997) showed the bight was sheltered from strong wind forcing throughout 
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the year. Wind stresses calculated from the SMB buoy and SCB buoy were relatively 
weak in comparison to the TB buoy, because SMB and SCB are within the bight and TB 
is on the western edge of the bight in the open ocean. Wind stress at all three buoys was 
strongest in April which was consistent with previous wind stress studies in the Southern 
California Bight (Winant and Dorman, 1997, Hickey, 1998, and Bray et al., 1999). 
During summer 2006 to summer 2007, the summer and fall wind stress at SMB and SCB 
were weaker than in the winter and spring. Winant and Dorman (1997) found that the 
wind stress tapers off in the summer and fluctuates in the winter due to propagating 


cyclones and anticyclones in the bight. 


The wind stress curl was compared to the alongshore currents to see if curl was 
the forcing the current reversals observed throughout the time-series. The cross 
correlation between the wind stress curl and the alongshore currents indicated that there 
was no statistical relationship between the current reversals and wind stress curl. In Bray 
et al. (1999), the correspondence between the current transport and either wind stress or 
wind stress curl was found to be not particularly good, except in the spring when the 
equatorward anomalies in transport near the coast were associated with equatorward wind 
stress anomalies. The coastal upwelling event observed at Eel Point was unrelated to the 


wind stress or the wind stress curl obtained from regional NOAA buoys. 
F. SIMPLE SLAB MODEL 


Near inertial/diurnal energy was observed in the ocean currents throughout the 
July 2006 to July 2007 time period. Since the frequency resolution of the wavelet spectra 
was too coarse to distinguish between diurnal or inertial currents, it was not possible to 
discern when the inertial motion actually occurred. Demodulated inertial/diurnal currents 
at 34 m peaked on April 4, 2007. However, the near inertial currents were not predicted 
at this time using a slab model driven by SCB buoy wind stress. The model was run with 
the demodulated currents from 38 m and showed that near inertial energy was confined to 
the near surface (shallower than 38 m). So, although it was possible to have near inertial 
frequency oscillations at the surface that did not penetrate down through the mixed layer, 


this could not be verified using the Eel Point observations because the shallowest 
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observed currents were below the speculated mixed layer. Note too that the winds used 


in the model were from a regional NOAA buoy and were not local winds. 


yZ 


Table 5. | Alonghore (AS) and cross-shore (CS) component of velocity for the currents off 
Del Mar, California and Eel Point, SCI. AS is positive toward 335° T (325 T) 
and CS positive toward 65° T (55 T). 


Standard Standard Principal axis 
Depth AS CS Deviation - AS | Deviation - CS orientation 
Season (m) (cm/s) | (cm/s) (cm/s) (cm/s) (deg T) 
[Summer] 60 | 13 | os | 6 | os | 335 
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Fall 
Winter 
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V. CONCLUSION AND FUTURE RECOMMENDATIONS 


The kinetic energy of the poleward flow along the coast off Eel Point was 
dominated by the clockwise rotating diurnal and semi-diurnal tidal variability. There 
were current flow reversals of low frequency flow observed each season, which were 
unrelated to regional wind stress and wind stress curl. There were also 2-3 and 3-4 day 
patterns of oscillations observed in the current velocity and stratification that were 
unrelated to regional wind stress and wind stress curl. Rotary spectra also revealed that 
there were dominant peaks of inertial energy at the near surface and near bottom depths. 
Wavelet power spectra were used to decompose the baroclinic tidal currents and revealed 
there were dominant modes of energy at the semi-diurnal and near inertial/diurnal period 
and that the strongest energy observed was at near the bottom. 

In April 2007, an oscillatory event was observed in the alongshore currents and 
temperature and salinity time-series. It was observed at the mid-depth downward through 
to near bottom depths in the alongshore velocity field. The disturbance had a dominant 
poleward flow with maximum velocity (30 cm/s) observed at mid-depth and a vertical 
phase velocity of 0.12 cm/s. These oscillations were rectilinear at mid-water and became 
counterclockwise with depth. It was speculated that this event may have been a coastally 


trapped wave. 


The wind stress within the bight was relatively weak throughout the time-series. 
Comparing the wind stress power spectra with along and cross-shore power spectra of the 
ocean currents, at the longer 6-10 day period observed in mid-November, there appeared 
to be some overlapping energy in both spectra. However, there was very little or no 
observed energy in the near inertial/diurnal and semi-diurnal or in the 2-3 day period in 
the wind stress power spectra that could be linked to the energy in the alongshore and 
cross-shore power spectra. There was no correlation between the wind stress curl and the 
alongshore currents. In conclusion, the wind stress and wind stress curl used in this study 
were not statistically correlated with the currents. Since there were no continuous local 
wind data available from SCI, the possibility of local (vs. regional) wind forcing could 


not be eliminated. 
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A. NAVAL RELEVANCE 


Current and future naval operations are focused on coastal/littoral regions around 
the world and often depend upon knowledge of local ocean conditions. Surprisingly, 
there is little known about the currents and circulation patterns around SCI. SCI/SCORE 
is owned and operated by Navy and supports training for the largest concentration of 
naval forces in the world. SCI/SCORE creates a realistic environment for tactical 
readiness training for military units and personnel as they deploy overseas. This thesis 
provides information on ocean currents and stratification on the continental shelf at SCI 
for SCORE and military participants to use with military exercises, regional 
climatologies, ocean engineering design studies, and marine mammal studies. Ocean 
current information is needed in all amphibious assault operations, mine hunting 
exercises, search and rescue (SAR), and special operations. The current and stratification 
descriptions are useful for all ASW operations. It is crucial when planning any operation, 
whether training or for war to have the environmental information for the operating area. 


This ocean data from Eel Point are meant to provide such data to operating naval forces. 
B. FUTURE RECOMMENDATIONS 


Obtaining local wind data closer to or at the mooring would help explain or 
establish the role of local wind forcing. Deploying another mooring at the southern end 
of SCI would allow the alongshore pressure gradient to be estimated. These 
measurements of pressure gradients could help explain the current accelerations and 
classify the physical character of observed waves. Since these observations will continue 
to July 2008, it may be possible that another year of current and temperature/salinity 
time-series can be used to see if the July 2006 to July 2007 period described here was 
typical. It would be beneficial to compare the currents to those produced by a local 
numerical ocean model such as ROMs Nowcast: _ http://ourocean.pl.nasa.gov/SCB. 
Finally, it would be useful to compare the variability of physical conditions observed at 


Eel Point with marine mammal vocalizations recorded near the moorings. 
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